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We calculated the Madelung energies of both the ground state and excited states in tetrathiafulva- 
lene - p-chloranil (TTF-CA) by taking into account intramolecular charge distribution. The distribu- 
tion is found to be significant in the neutral-ionic (NI) transition. In the ionic phase, the Madelung 
energy depends more strongly on the intermolecular distance perpendicular to the 7r-stacking chains 
than on that along the chains. This indicates that simple single-chain models neglecting interchain 
electrostatic coupling are not adequate. The gain of the Madelung energy due to dimerization is 
concluded to be small compared with the other structural changes. We also calculated the formation 
energy of excited state domains, which appear in the primary process of the phase transition. A 
one-dimensional excited domain has the smallest energy among the possible domains with the same 
number of molecules when the domain is small, and the energy per molecule is considerably reduced 
in increasing the domain size. It is consistent with the experimental suggestion that a large number 
of excitations were generated by only one photon. 
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I. INTRODUCTION 

Controlling properties of materials by various external 
stimuli has been attracting attention from the viewpoint 
of fundamental researches and applications. An organic 
compound tetrathiafulvalene - p-chloranil, TTF-CA, is 
one of the typical model materials ||^-@| , where a phase 
transition accompanied by the change of dielectric re- 
sponse |6|] is induced by applying pressure, or varying 
temperature, or irradiation of light. Torrance et al. first 
discovered the pressure-induced neutral-ionic (NI) phase 
transition of TTF-CA [n], followed by the demonstration 
of the temperature- induced one g] . At high temperature 
(or at low pressure), both TTF and CA molecules are 
nearly neutral in their electronic charges (N-phase). On 
the other hand, at low temperature (or at high pressure), 
they are quasi-ionic (Tphase) due to electron transfer 
from TTF to CA. The degrees of charge transfer p were 
estimated to be about 0.3 in the N-phase and 0.7 in the 
Tphase 0. The photo- induced changes relevant to the 
NI transition were also reported 0,||,^. The I-phase is 
partially converted to the N-phase by irradiation of laser 
pulse with photon energy 2.0-2.5 eV. It was reported that 
about 300 pairs of TTF-CA molecules were converted by 
one photon |^ , suggesting that this phenomenon involves 
a cooperative effect with large nonlinearity. 

The planer TTF and CA molecules shown in Fig. Ill 
form mixed-stacking chains along the a-axis in the crys- 
tal PI . The temperature-induced phase transition as well 
as the pressure-induced one is understood in terms of 
changes in electrostatic Madelung energy by lattice ex- 
pansion or contraction. A dimerized structure was ob- 
served in the I-phase ||] . The origin of the dimerization 



was also ascribed to the electrostatic interaction in pre- 
vious theoretical works p^|-p^ . 

In most of the previous works, TTF-CA was treated as 
a quasi one-dimensional chain along the mixed-stacking 
a-axis with weak or no interchain coupling with respect 
to the Madelung interaction as well as to the overlapping 
of the TT-orbitals P,0,H. That is, the electrostatic inter- 
action between chains was assumed to be much smaller 
than that within the chain Pq]. This assumption was 
justified in the approximation where each molecule was 
replaced by a point charge at its center of mass. Here- 
after this is called the point molecule approximation. The 
nearest neighbor molecules of TTF or CA are the same 
along the 6-axis and the different along the b-\-c direction 
|l6| . As a result, a repulsive interaction along the 6- axis 
and an attractive one along the b -\- c direction are ex- 
pected to cancel out each other iQ. However, this idea 
is not consistent with the change observed in the lattice 
constant [^ . The lattice constant b becomes small at the 
NI transition with lowering temperature. This implies 
an attractive interaction along the &-axis. In contrast, a 
change in the lattice constant a is very small, although a 
magnitude of the interaction changes at the transition. 

In the present paper, we will demonstrate the inade- 
quacy of the point molecule approximation. In the calcu- 
lation of Madelung energy, we take account of the charge 
distribution on atoms inside each molecule derived with 
an ab initio quantum chemical method flT]] . We call this 
approach the point atom approximation. Another pur- 
pose of this paper is the elucidation of the primary pro- 
cess of the phase transition. We estimate the formation 
energy of an excited-state domain, e.g. a neutral do- 
main in the ionic phase. This will provide information 



about the growing process of the excited domain and 
inter-domain interactions. 

This paper is organized as follows: The charge dis- 
tribution of each molecule calculated with an ab initio 
technique is presented in Section II. In section III, the 
Madelung energy of the ionic phase is calculated with 
various modifications of the crystal structure to show the 
inadequacy of the point molecule approximation. Section 
IV deals with the formation energy of the excited-state 
domain. Finally, we give concluding remarks in Section 
V. 



II. CHARGE DISTRIBUTION INSIDE 
MOLECULES 

We calculated the intramolecular charge distribution 
in each molecule using the ab initio quantum chemical 
method [H. The calculations were performed both for 
neutral and monovalent states. The molecular structures 
were extracted from the crystal structure of TTF-CA 
determined at 90 K, in which the molecules have the 
inversion symmetry B. The restricted and unrestricted 
Hartree Fock methods were respectively used for the neu- 
tral and monovalent states with the 6-31G* basis function 
set. The charge density on each atom was calculated with 
the MuUiken population technique. The obtained den- 
sity is consistent with the previous studies |jl^,^. The 
charge distribution at an observed fractional charge p is 
evaluated by a linear interpolation between the complete 
neutral and monovalent states. 

The results are shown in Fig. H. The large electronic 
polarization is found on both the molecules. For exam- 
ple, the carbon atom next to an oxygen atom in the 
CA molecule has a large positive charge even if the to- 
tal charge of the molecule is negative. The derived in- 
tramolecular charge distribution is used for the point 
atom approximation in the following sections. 



III. STRUCTURE DEPENDENCE OF 
MADELUNG ENERGIES 

In this section, we report the Madelung energy, Eyi, 
of the ionic phase of TTF-CA using the Ewald method. 
We examine how the electrostatic energy changes with 
the following deformations of the crystal structure: the 
variation of the lattice constants, the molecular rotation 
around its center of mass, and the dimerization along 
the a-axis. The intramolecular atomic coordinates and 
the charge distribution are fixed because intramolecular 
deformations are smaller than intermolecular ones. This 
means that only the spacing among molecules is varied. 

Figure shows the variation of Em against the uniax- 
ial deformation along the a and b axes from the crystal 
structure observed at 40 K ||^ . The results are quite dif- 
ferent between the two approximations. The lowering of 



Em due to the contraction of the lattice constant a is 
suppressed in the point atom approximation. As for the 
case of the 5-axis, the slope of Em against lattice con- 
stant is different in the two approximations. In the point 
molecule approximation, the slope is slightly negative, 
which agrees with the previous proposal of the repulsive 
electrostatic interaction along the 6-axis |lq,Qq]. In con- 
trast, in the point atom approximation, Em is decreases 
in contracting along the 6-axis. In this approximation, 
the attractive interaction is stronger along the 6-axis than 
along the a-axis, which is consistent with the observation 
that the lattice constant along the 5-axis is much short- 
ened at the phase transition temperature in going from 
the N-phase to the I-phase ||] . With varying the lattice 
constant along the c-axis, we found a tendency similar to 
the case of the 6-axis. 

Figure y shows the electrostatic energies between 
neighboring molecules in the a6-plane. The attractive 
energy between TTF and CA aligned along a/2 -I- 6, 
defined as 6'-axis, is found to overcome a repulsive one 
between CA's along the 6-axis, while the intermolecular 
distance measured in terms of the center of mass of each 
molecule is larger in the former pair than in the latter. 
This indicates the inapplicability of the point molecule 
approximation. The discrepancy is mainly caused by the 
hydrogen bond along the 6'-axis. One of the hydrogen 
atoms of TTF is very close to one of the oxygen atoms 
of CA. This hydrogen bond makes a large contribution 
to the interchain attractive interaction. 

The inclination of the molecular planes against the 
stacking axis also has some influences on Em- At the 
NI transition, the molecules rotate mainly within the ah- 
plane. We define the inclination angle (j) as shown in the 
inset of Fig. |[ The angle (j) of the CA molecule (0ca) 
in the I-phase (40 K) is found to be larger by 1.4 degree 
than that in the N-phase (90 K), and found that of the 
TTF molecule ((/ittf) in the I-phase is also larger by 0.7 
degree. These rotations alter the effective intermolecular 
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distance R cos in the inset of Fig. 4 . 
Em on i? cos 0c A is shown in Fig. 4. 
angle of TTF is set to be half that of CA to match the 
observed changes of 0cA and 0ttf- The other geomet- 
rical parameters are fixed to the crystal structure deter- 
mined at 40 K. Em is an increasing function of i? cos 0cA- 
This means that the Madelung energy is stabilized by the 
decrease of the effective intermolecular distance. The 
energy difference due to the molecular rotation at the 
NI transition is 80 K/formula, which is nearly the same 
magnitude as the effect of the contraction of the lattice 
constant a by 1% as shown in Fig. y. This implies that 
the enlargement of an attractive interaction along the a- 
axis in the I-phase makes the effective distance between 
molecules to be short instead of the contraction of the 
lattice constant. 

The dimerization observed in the ionic phase has often 
been considered as an essential feature of the NI transi- 
tion [ pTp^ . We also calculated the stabihzation energy 
due to the dimerization. The gain of the Madelung en- 



ergy is 2.3 K/formula for a displacement of CA by 1 % 
along the a-axis, which is comparable to the observed dis- 
placement in the ionic phase. The stabilization energy is 
much smaller than the energy changes due to the other 
deformations studied in this paper, as well as than the 
observed phase transition temperature, 81 K. 



IV. FORMATION ENERGY OF THE EXCITED 
DOMAIN 

In this section, we consider an excited-state domain 
in the ground state, e.g. a neutral domain in the ionic 
phase, in connection with the photoinduced NI tran- 
sition. We assume that a pair of excited TTF-CA 
molecules along the a-axis is a constituent unit of the 
domain, because the transfer of a tt electron is likely to 
occur between TTF and CA along the stacking axis. The 
constituent unit and some examples of the domain are 
given in Fig. ^(a). We calculated the formation energy 
^EX p£ -j-j^g domain in an extended scheme of the point 
atom approximation. Here we show the result with the 
crystal structure at 90 K for simplicity and for avoiding 
the influence of the dimerization. E^^ is divided into 
two parts, 
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where A^ex represents the number of the constituent units 
in the domain. E^ originates from the Madelung energy 
described as follows: 
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where GS and EX indicate the ground state and the ex- 
cited state of the constituent unit, respectively. E^ de- 
notes the Madelung energy in the excited domain, and 
E"j^^ is that between the domain and the surrounding 
ground state. Group D represents the set of atoms in the 
excited domain, and i and j specify atoms, q" represents 
the charge of the i-th atom in the a-state, and Rij means 
the distance between the i-th and j-th atoms. The sec- 
ond term of E^^ in Eq.yJ represents all the contributions 
other than the Madelung energy, for example, electron 
affinity, ionization potential, and electron transfer ener- 
gies. AEu is defined as E^^ — E^^, where E^ represents 
those contributions in the a state per constituent unit. 
Instead of calculating E" directly, we assumed that the 
N- and I-phases are degenerate in energy. In this case, 
AE'u is equal to the difference of the Madelung energies 
per unit between the two phases. In the expression of 
E^, we find with a straightforward calculation that the 
formation energy does not change by interchanging gp^ 



and qf^. This means that the N-domain in the I-phase 
has the same energy as the I-domain of the same shape 
in the N-phase. 

We investigated E^^ of the one-dimensional (ID) do- 
mains with various shapes, specifically along the a, b, b' , 
and c-axes with the number of pairs ranging from 1 to 
50. Figure H(b) shows E^-^ of ID domains along the a 
and 6'-axes, as a function of the domain size. For com- 
parison, the results using the point molecule approxima- 
tion are also shown. The formation energy is found to 
be sensitive to the direction of the domain growth. It 
is the lowest in the direction along the a-axis among all 
the types; that for the 50 excited pairs along the a-axis 
is at 0.016 eV per pair, which is only 1/8 as large as 
that of a single excited pair. Such a drastic reduction of 
the formation energy accords with the experimental sug- 
gestion that a large number of excited pairs are created 
by a single photon in the photoinduced change p]. The 
formation energy for the ID domain along the 6'-axis is 
much reduced by taking account of the molecular struc- 
ture as for the comparison between the two approxima- 
tions. This indicates, again, that the interchain coupling 
is much stronger in the real situation than in the point 
molecule approximation. 

E^^ of the two- or three-dimensional excited domains 
were also calculated. In Fig. ^, we show the results for 
two-dimensional domains within the (a, b') plane. The 
excited domain has a shape of parallelogram consisting 
of A'ox excited pairs spanned by Na and A^;,' units along 
the a and 6'-axes. It is found that the excited domain 
with Nbi = 1, i.e. the ID domain along the a-axis, is the 
most stable among the domains of the same size in the 
case of A'ex=12 or 18. This remains true as long as the 
number of excited pairs is less than twenty. This suggests 
that excited domains easily grow along the a-axis in the 
primary stage of the phase transition. 



V. CONCLUDING REMARKS 

Here we comment on the electron transfer interaction 
that has been neglected in our calculations. Since the 
transfer integral of tt electrons between chains is consid- 
ered quite small, the essential role of the attractive elec- 
trostatic interaction between chains would be unchanged 
even if we take account of the itinerancy of electrons. 
On the other hand, the electron transfer is considered 
to have some infiuence on the electronic structure within 
chains. The fractional degree of the charge transfer that 
we assumed in this paper is ascribed to this interaction. 

The growth process of the excited domains in chains 
will also be affected. As well as the electrostatic inter- 
action, the electron transfer is expected to accelerate the 
domain growth along the mixed-stacking chain. With the 
intermolecular electron transfer energy, the dynamics of 
the photoexcited state has been investigated by means 
of a one-dimensional modified Hubbard model ||20|l, in 



which the electrostatic interaction between molecules was 
neglected. For more quantitative description of the dy- 
namics, it is necessary to develop a theory that includes 
both the electrostatic interaction with the intramolecular 
charge distribution and the electron transfer interaction. 

Next, we discuss the driving force of the photoinduced 
phase transition. We have calculated the effect of the lat- 
tice constant expansion or contraction on the Madelung 
energies. The 0.7 % contraction of the lattice constant 
b, which was the observed difference between 40 K and 
90 K, lowers Em about 20 meV per TTF-CA pair. The 
energy gains due to the molecular rotation within the 
a6-plane and the dimerization along the a-axis is about 
7 meV and 0.2 meV, respectively. On the other hand, 
the formation energies E^^, derived without any lattice 
relaxations, of an excited TTF-CA pair and the excited 
domain of 50 pairs are respectively 135 meV and 16 meV 
per pair. Such a drastic reduction of the formation energy 
is caused by the electrostatic interaction. The change in 
the Madelung energies is likely to be essential in the pri- 
mary process after the photoexcitation. The change of 
the lattice constant b also becomes important in increas- 
ing the excited domain size, where E^^ becomes com- 
parable to the lattice relaxation energy. The role of the 
lattice constant variation was discussed for the coopera- 
tive photoinduced phase transition of spin-crossover com- 
plexes pi[-p3|. Similar discussion may also be required 
in TTF-CA. 

In summary, we calculated the electrostatic energies in 
both the ground and excited states of TTF-CA. It has 
been concluded that the point molecule approximation, 
where each molecule is replaced by a point charge at its 
center of mass, is inadequate for this system. For the 
ground state, we have demonstrated the crucial effect 
of interchain interactions, especially between TTF and 
CA molecules aligned in the a/2 -I- 5 direction. For the 
excited states, we found large nonlinearity in the size 
dependence of the formation energy of the excited-state 
domain. The growth of the domain will proceed along the 
a-axis in the primary process after the photoexcitation. 
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FIG. 1. Structure of TTF and CA molecules. The cal- 
culated charge distributions on each atom are also shown for 
p=0.3, and for p — 0.7 in parentheses. 
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FIG. 2. Shifts of _Em due to the lattice constant variation. 
The abscissa axis represents each lattice constant normalized 
by that observed at 40 K, £4ok[9]. The results using the point 
atom approximation and with the point molecule approxima- 
tion are represented by 'atom' and 'mol.', respectively (see 
text). 
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FIG. 4. Shift of _Em due to the molecular rotation from 
the structure at the I-phase (40 K), where R and 0ca repre- 
sent the distance between the center of mass of neighboring 
CA molecules along the a-axis, and the molecular inclination 
angle for CA molecules defined in the inset, respectively. The 
angle 'F and 'N' indicate the inclination angles (/jca observed 
at the Fphase (40 K) and at the N-phase (90 K)[9]. 
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FIG. 3. Molecular alignment within the afe-plane. The in- 
termolecular electrostatic energy is shown for each pair. The 
distances between their centers of mass are also indicated 
along in parentheses. 



(a) 



b'=a/2+b 




(b) 



0.25 



^0.10 




■ atom:fl-axis 
:^'-axis 
- mol. :fl-axis 
:4''-axis 



10 20 30 40 
N(pair) 



50 



0.045 

0.040 

|h 0.035 h 
> 

^0.030 

I 0.025 

0.020 

0.015 



I -r— 



N=\ —m—N =12 

D ex 

N=2 —A—N =18 

o ex 

N=3 —m—N =30 




10 20 30 

A^fpair) 



40 



FIG. 6. Formation energy of two-dimensional excited do- 
main in the (a, b') plane, where Na, Ny , and Ncx{= Na x Ny) 
represents the number of TTF-CA pairs along the a-axis in 
the domain, that along b' — a/2 + b, and total number of pairs 
in the domain. Domain sizes Na x Ny explicitly indicate for 
iVcx=30. 



FIG. 5. (a) Examples of the structure of ID domains. 
A constituent unit of the excited domain is represented by 
'Unit', (b) The formation energy of the ID domains along 
each axis. The results using the point atom approximation 
and the point molecule approximation are represented by 
'atom' and 'mol.', respectively (see text). 



